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Abstract 

Let Q\ (t, y) be the number of people present at time t with y units of remaining service time 
in an infinite server system with arrival rate equal to A > 0. In the presence of a non-lattice 
renewal arrival process and assuming that the service times have a continuous distribution, 
we obtain a large deviations principle for Q\ (•) /A under the topology of uniform convergence 
on [0, T] x [0,oo). We illustrate our results by obtaining the most likely path, represented as 
a surface, to ruin in life insurance portfolios, and also we obtain the most likely surfaces to 
overflow in the setting of loss queues. 

1 Introduction 

The asymptotic analysis of queueing systems with many servers in heavy-traffic has received sub- 
stantial attention, especially in recent years. Among the earliest references that come to mind 
in connection to this topic is the work of [7] on heavy-traffic limits for the infinite-server queue. 
Another highly influential paper in the area is [6] in the context of many server Markovian queues, 
which introduced a scaling that is now known as the "Quality and Efficiency Driven" regime. The 
ideas in these papers have fueled more recent results in the asymptotic analysis of many server 
systems such as: [12], [E], [13], |D], [ID], in the setting of many server queues, and [5], [2], [JT], [Ti] , 
in the setting of the infinite server queue. The asymptotic analysis of queueing systems with many 
servers has been motivated by applications in service engineering, in particular in the context of 
call centers and health-care operations. Another set of application areas that is also very relevant, 
but that is infrequently mentioned in the analysis of many server systems is that of insurance math- 
ematics. It is clear, for instance, that a portfolio of insurance policies can be directly modeled as 
an infinite server system; casting insurance portfolios in this framework is particularly appealing in 
the setting of life insurance as we shall illustrate in Section [5] 

So far most of the asymptotic analysis of many server systems has concentrated mainly on 
fluid and heavy-traffic approximations. Meanwhile, the literature on large deviations analysis for 
many server queues is not as extensive as that of fluid or diffusion approximations; despite the fact 
that it is clearly of interest to understand the large deviations behavior of these types of systems. 
For instance, consider the consequences of dropping calls in an emergency call center or being 
unable to satisfy the demand for critically ill patients in the context of health-care applications. 
In the insurance setting, for instance, it is of interest to estimate ruin probabilities and, perhaps 
even more importantly, understanding the most likely path (or set of paths) to ruin. Risk theory 
typically concentrates on ruin probabilities for aggregated models, such as the classical ruin model 
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(see PQ); the results in this paper, as we shall illustrate, provide a systematic way for assessing ruin 
probabilities for a class of bottom-up models. 

Our main contribution in this paper is to provide the first sample-path large deviations analysis 
of the state descriptor of the infinite server queueing model in heavy-traffic (i.e. as the arrival 
rate increases to infinity without introducing any scaling on the service times). The statement of 
our main result, which is given in Theorem [TJ features a convenient representation of a good large 
deviations rate function, under a strong topology, that later we use in some applied settings of 
interest. For instance, we compute the most likely path to overflow in a loss system, and also the 
most likely path to ruin for a life insurance portfolio that embeds an infinite server queue with a 
particular service cost structure. It is important to emphasize that our results take advantage of a 
convenient representation of the system's description that facilitates the representation of the rate 
function; detailed discussion on this system's representation is given in Section 2.1 Previous large 
deviations analysis of the infinite server queue has concentrated on queue-length characteristics 
only; see, for instance, [lj who develops large deviations for marginal quantities in the case of 
renewal arrivals, and [15] who develops sample path large deviations for the queue length process 
of infinite server queues in tandem in the case of Poisson arrivals. 

Our large deviations analysis complement results on fluid analysis and diffusion approximations 
recently obtained for infinite server systems. For example, [TT] have shown that the state descriptor 
of the infinite server queue, suitably parameterized in terms of a two-parameter stochastic process, 
converges after centering and re-scaling to a Gaussian process; see also |14| who interpret the state 
descriptor of the infinite server queue as a measure valued process acting on the space of tempered 
distributions. These recent results, in turn, extend prior work by [5] in the context of discrete and 
bounded service time distributions, and [2] for the case of Poisson arrivals. 

The analysis of the infinite server queue is important as it serves as a building block for other 
models of interest. For instance, in the setting of loss models one can clearly couple the loss 
systems with associated infinite server systems, and in the setting of many server queues |13j shows 
how one can precisely understand queues with multiple servers as a perturbation of infinite server 
queues. Furthermore, the infinite server model is a classical model in queueing theory that serves 
as direct model in important applications. Of particular interest to us, as mentioned earlier, are 
the applications to insurance mathematics. 

The rest of the paper is organized as follows. In Section [2] we introduce our problem setting 
and provide a statement of our main result. This is a fundamental section and it is divided into 
three parts. We first shall introduce our assumptions and define our notation. Then we provide 
the precise mathematical statement of our result, and, finally, we will provide a heuristic argument 
that allows us to gain some intuition behind our main result. The next two sections then provide 
the proof of our main result. We first show our result for bounded service times in Section [3] Then, 
in Section |4j we apply a truncation argument to extend our result to unbounded service times. 
Finally, in Section [5] we apply our result to computing the most likely paths to rare events in the 
setting of loss queueing systems and also in the setting of ruin probabilities for large life insurance 
portfolios. 



2 Assumptions, Notation, and Main Result 

The purpose of this section is threefold. First we shall clearly state our assumptions and introduce 
necessary notation for our development. Second, we shall explain the main large deviations result 
and provide a heuristic derivation of the rate function that we obtain. Finally, we shall provide a 
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road map for the strategy behind the proof which will be presented in subsequent sections. 
2.1 Assumptions and Notation 

We shall describe an underlying system corresponding to an arrival rate A. We call the system with 
A = 1, i.e. one customer per unit time, our "base system"; eventually we shall send A to infinity in 
our asymptotic analysis. We collect our assumptions as follows. 

Assumptions and notation concerning the arrival process. For the base system, we assume 
the interarrival times are non-lattice, i.i.d., positive random variables (U n : n > 1) with finite 
exponential moments in a neighborhood of the origin; in precise words, k(6) := logEe eUn < oo for 
some 9 > 0. In our A-scaled system, the arrivals come A times faster (i.e. the n-th interarrival 
times becomes U n /X). The associated logarithmic moment generating of the A-scaled service times 
is then K\(Q) := \ogEe eUn / x = k{6/\) < oo for some 9 > 0. 

The time at which the n-th arrival occurs in the base system is A n = U\ + ... + U n for n > 1. 
We simply define Aq := and then let ./V (t) := maxjra > : A n < t} be the number of arrivals that 
have occurred up to time t in the base system. It is important to keep in mind that ./V (•) increases 
by one unit at discontinuity points since we are assuming that the J7 n 's are positive. 

Eventually, we shall increase the arrival rate, so it is sensible to define N\ (t) := N (\t). 

Define the so-called infinitesimal logarithmic moment generating function for the arrival process 
via iPn{9) = —k~ 1 (—9) (see [1]). This definition is motivated by the fact that 

lim A" 1 log£exp (9[N X (t + 6)- N x (i)]) = %f> N (9) 5 

t— >OD 

for any 5 > 0. Since the U n 's are positive with probability one we have that ipM (•) is continuous 
and strictly convex on the positive line. We also assume that ipN (•) is continuously differentiable 
throughout R. This assumption is satisfied for most arrival processes, certainly for interarrival 
times that are strictly positive and such that sup{ft (0) : k(9) < oo} = oo. 

Assumptions and notation concerning the service times. We assume that the n-th customer 
that arrives to the base system (i.e. at time A n ) brings up a service requirement of size V n . The 
sequence (V n : n > 1) is assumed to be i.i.d. We write F(x) = P(V n < x) to denote the associated 
distribution function evaluated at x, and set F(x) := 1 — F(x) to be the tail distribution. Moreover, 
we assume that F (■) is continuous. 

Two-parameter representation of system status. Let Q\ (t, y) denote the number of cus- 
tomers who arrived before or at time t and leave after time y in the A-scaled system. In other 
words, 

Qx {t, y) = Qx (0, y - t) + 1 (Yn + ^n/A > y) . 

71=1 

We shall assume that the system is initially empty at the beginning. This is done for simplicity. 
Since we have infinitely many servers, we can incorporate the initial configuration by keeping track 
of its evolution independently of what occurs subsequently. Given our assumption of an initial 
empty system we then have that Q\ (0, u) = for all u > 0. Note that for t G [0, T] and y > 0, 



Qx (t,y)=Q x (tAy,y)+N x {t)-N x (tAy) 



(1) 
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It is worth comparing the current system representation with the more common one involving 
the quantity Q\ (t, u) defined as the number of customers in the system currently at time t who 
have residual service time larger than u > 0; more precisely, Q\(t,u) = Q\(t,u + t). These two 
system representations are equivalent in the sense that (Q\ (t,u) : t G [0,T],u > 0) encodes the 
evolution of the infinite server systems and thus, such evolution can be used in principle to retrieve 
(Qx(t,u) : t G [0,T],u > 0). We have chosen the representation based on Q\ to facilitate the 
representation of the rate function; a more detailed discussion is given towards the end of Section 



2.2.1 In addition, the representation based on Q\ allows to obtain a rich large deviations principle 
to which one can apply the contraction principle directly to several continuous functions of interest. 
For instance, it follows immediately that the arrival process N\ (t) = Q\ (t, 0), and the departure 
process, D\ (t) := N\ (t) — Q\ (t,t) are continuous functions under the topology that we consider 
(and that we shall discuss in the next paragraphs) . More applications of the contraction principle 
will be discussed in Section [U 

Discussion about the topological space. Let V = {(t,y) : < t < T, y > 0} and let us write 
|| • ||c to denote the supremum norm over any set C. The space of functions that we consider for our 
large deviations principle shall be denoted by L +>0O (T>) and it corresponds to bounded functions 
with domain in D, such that x (0,u) = for u > 0, x (t, •) is non increasing, and x (t, ■) vanishes 
at infinity. We will develop the large deviations principle for the family of stochastic processes 
(<3a/A : A > 0) on the space L +)OQ (T>) endowed with the topology generated by the supremum 
norm. Following [3j p. 4, the probability measures in path space in our development are assumed 
to have been completed. 

Our large deviations principle for Q\/\ immediately implies in particular a large deviations prin- 
ciple in the Skorokhod topology in the space D D i 0jOO ^[0,T] which is the space of right-continuous- 
with-left-limits (RCLL) functions x, with domain on [0, T], that take values on the space of RCLL 
functions taking values on M. That is, on each time point t in x = (x (t) : t £ [0, T]) G -D_d r [o,oo) [0, T] 
is a function x (t) G -Dr[0, oo). This is precisely the topology considered in who also provide a 
discussion on the benefits of using this topology relative to other natural (but weaker) alternative 
options (see Section 2.3 in [TT]). 

An alternative approach that one might consider given the available results on functional weak 
convergence analysis of the infinite server queue, such as |14j . is to interpret the space descriptor 
of the infinite server queue as acting on the space of tempered distributions. We believe, however, 
that this approach, although elegant, has important limitations in terms of assumptions and the 
class of functions to which the contraction principle can be directly applied to obtain other large 
deviations principles of interest. 



2.2 Statement of Our Main result. 



We are now ready to state our main result. Let q :- 
q G AC + (V) if the following conditions hold: 
i) q is absolutely continuous on V, and 



d 2 



dt dy 



q(t,y) 



: (t,y) : (t,y) G V) G L +j00 (V). We say that 



dydt < oo, 



ii) d 2 q(t,y)/dt dy = almost everywhere for (t,y) G {(t, y) : < y <t < T}, 

iii) q(0,y) = for y > 0. 

If q G AC + (V) , then we let I (q) be defined via 
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sup 

9( v )eC(,[0,T]x[0,oo) Jo 



9(t,y-t) 



d 2 
dt dy' 



(t, y))dy-il) N [ lo, 



'■' y) dF(y) 



dt (2) 



where C7j[0, T] x [0,oo) is the set of all bounded continuous functions on [0, T] x [0,oo). On the 
other hand, if q £ L +j00 (P) fails to satisfy any of the conditions i) to hi), simply let / (q) = oo. 
We now can state our main result. 



Theorem 1. Under the set of assumptions discussed in Section 2.1, (Q\/X : A > 0) satisfies a large 
deviations principle with good rate function I (•) on the space (L +j00 (P), ||-||-p). In precise terms, 
for each open set O we have that 



!imA->oo lQ g t p (<2a/A G O) > - inf I (q) 
A qeu 



and for each closed set C 



lim^so log -P(Qx/X E C) < - inf / (q) . 

A q&C 



As an immediate corollary we obtain, as mentioned earlier, a large deviations principle for 
(Qa/A : A > 0) under the Skorokhod topology in the space D D r 0jOO ^[0,T], discussed in the previous 
section and introduced in 

We shall explain the strategy behind the proof of Theorem [TJ First, w e sha ll introduce an 
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We will show 



auxiliary continuous process, Q\/X, which shall be defined later in Section 
that Q\/X is exponentially equivalent to Qa/A. Second, in addition to the assumptions imposed 
in Section |2.1| we will assume that there exists a deterministic constant M £ (0, oo) such that 
P (V n £ [0, M]) = 1. In the third and last part of the argument we will relax this truncation 
assumption. 

In turn, the first part of the argument (i.e. assuming truncation) is divided into several steps. 

The first step consists in developing the large deviations principle for Q\/X with rate / (•) 
under the topology of pointwise convergence using the Dawson-Gartner projective limit theorem. 
The second step involves showing that Q\/X is exponentially tight as A — > oo under the uniform 
topology on the compact set [0,T] x [0, M\. The third and last step involves lifting the large 
deviations principle to the uniform topology. 

During the second part we introduce an approximation scheme that proceeds by ignoring the 
customers that arrive to the system with a service time larger than K. Using a coupling argument, 
the process that is obtained using this scheme is shown to be a good approximation to the original 
system for the purpose of computing large deviations probabilities. 

However, before we do this let us provide a heuristic argument in order to guess the form of the 
rate function. Later we will explain what are the technical difficulties that need to be addressed. 



2.2.1 Guessing the Rate Function: A Heuristic Approach 

One can take advantage of the point process representation of the input process (i.e. the arrivals 
and the service times represented as a marked point process). Let us start with the case of Poisson 
arrivals. We shall briefly explain how to adapt the development that follows to the more general 
case of renewal arrivals. 
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Consider the scaled system with arrival rate A and suppose that F (•) has a density / (•). The 
amount of customers that arrive during the time interval [t,t-\-dt] that bring a service requirement 
of size [r, r + dr] is denoted by the quantity M\ (t + dt,r + dr), which is governed by a Poisson 
distribution with rate A/ (r) dtdr. It follows then by elementary considerations involving the Poisson 
distribution that M. \ (t + dt, r + dr) /A satisfies a large deviations principle in the real line. In 
particular, we formally obtain that 

P(M X (t + dt,r + dr) /A w fj, (t, r) dtdr) = exp (-AJ (fi (i, r)) dtdr) , 

with 

J O r )) = sup [77 (i, r) /i (t, r) - ip N {V (*> 0) / (?")]> 

r;(t,r) 

and V'iv (??) = exp (77) — 1. The supremum above is obtained formally with 77* (t, r) := log(/i [t, r) / f (r)). 

So, by pasting independent regions of the form [t, t + dt] x [r, r + dr] together one expects that 
the Poisson random measure A4\ (•) /A would satisfy a large deviations principle under a suitable 
topology, so that 

P ( M\ (Ax B) /\& fj, (t, r) dtdr, for a large class AxB a exp (-AJ (/i)) 



with 



J Ox) 



[0,T]x[0,oc) 



sup 



[77* (i, r) /i (t, r) - tpN (v* (t, r)) f (r)]dtdr 



t)(;-)eC b [0,T\ x [0,oo) J [0,T] x [0,oo) 
Now, observe that for all y > t 

Qx (t,y) 

' [0,t] x [y— s,oo) 

and (Qa (0, y) = for y > 0. Now, consider 



[77 (i, r) /j (i, r) - (t/ (*, r)) f (r)]dtdr. 



M\ (s + ds,r + <2r) 



Note that 



dydt 



t roo 



v (s, r) drds. 



(3) 



(4) 



(5) 



Jy-s 



q(t,y) = -v (t,y- 1) 



Therefore, if q(-,-) is absolutely continuous and q(0,y) = for y > 0, so that representation 
([5| is applicable, one can formally compute the rate function of Q\ (•, •) / A evaluated at q(-,-) by 
evaluating J (77) with 

d 2 

^ r) = ~dydt^ V) 
for s G [0,T] and r G [0, 00). In particular, this analysis yields that I (q) should satisfy 

[77 (s, r) /j (s, r) - Vat (77 (s, r)) / (r)]drds 



sup 

r7( v )6C 6 [0,T] x [0,oo) J [0,T] x [0,oo) 



sup 

r?( v )GC*,,[0,T] x [0,oo) J [0,T] x [0,oo) 



77 [s,r) 



J) 2 
<97/<9i 



g(s, s + r) - t/^at (77 (s, r)) / (r) 



sup 

r?( v )eC b [0,r] x [0,oo) J [0,T] x [s,oo) 



7/ (s, 77 — s) 



a 2 



<97/<9t 



drds 



q(s, u) - (exp(7/ (s, u - s)) - 1)/ (77 - s) 



duds, 
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which is, of course, equivalent to ^ in the Poisson case assuming the existence of a density / (•) 
for the distribution of the service times. The previous form of the rate function was heuristically 
obtained assuming that y > t. However, since all the information of the infinite server queue is 
contained in the evolution of (Q\ (t,y) : t G [0, T],u > 0) and Q\ (t,u) = Q\ (t,t + u), we must 
have that the rate function should be specified only over q(t,y), such that t G [0, T] and y > t so 
that d 2 q(t, y)/dydt = 0if0<y<t<T. 

For the non-Poisson case one can argue using renewal arguments. We need to compute the 
log-moment generating function of the vertical strip (t + dt, + dr) : 1 < i < n), where 
r i < r 2 < ■■■ < r n for an arbitrary partition (r^ : 1 < % < n). We obtain, using elementary properties 
of the multinomial distribution together with an application of the key renewal theorem as in [3], 



E 



E 



E 



N(\(t+dt))-N(Xt)' 



exp ^ 9 (i, n) M\ {t + dt, n + dr) 
\i=i / 

( exp (9 (t, n)) P {V x G [r h r t + dr}) 

\i=l / 

exp \[N (\{t + dt)) - N (Xt)] log ex P ( r i)) P (^i G r * + dr ]) 
exp ( A^at (log (^exp(0(t,r i ))P(y 1 G [r h n + dr]) ) J +o(A)J 



vi=l 



as A — > oo. 

So, by pasting together vertical strips (i.e. ranging the parameter t) we obtain that the family 
of random measures M\ (•) /A is expected to satisfy a large deviations principle under a suitable 
topology with rate function 



sup 

e(;-)eC b [0,T]x[0,oo)J[0,T] 



oo / poo 

0(t,r) fj,(t,r)dr-ip N I / exp (9 (t,r)) dF (r] 
o VVo 



The rest of the formal analysis proceeds similarly as in the Poisson case. 

The formal argument just outlined, even if heuristic, suggests a potential approach to developing 
sample path large deviations for Q\/X. Namely, first develop a large deviations for the random 
measures M\ (•) /A, and then apply the contraction principle to obtain the desired large deviations 
result for Q\/X. This approach, although intuitive, will not be followed in our development. We 
found it easier to directly work with the topology that we wish to impose. Part of the problem 
involved in making the argument based on random measures rigorous in the setting of the topology 
that is of interest to us is that indicator functions are not continuous, so the contraction principle is 
not directly applicable if one is to endow the space of measures with the weak convergence topology. 
Of course, one can proceed by trying a different topology (stronger than weak convergence) or by 
trying to use the extended contraction principle. However, the technical development, we believe, 
would end up being more involved than the direct approach that we will follow. 

An additional concern that might arise at this point is our selection of Q\/X in order to represent 
the system status; as opposed to Q\/X, which might appear more natural at first sight. Let us 
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explain why Q\/\ is a more convenient object to consider. Note that if q (s, r) = q(s, s + r), then 



c) 2 



9s dr 



'(s,r) 



& 1 



ds dr 

a 2 

ds dr 



8 2 

q(s, s + r) + 7r^g(s, s + r) 

a 2 

g(s,s + r) + ^2?(s,r) 



so 



9 2 



ds dr 



d 2 d 2 



Since Qa (£j u) /A = Qa (£, it + i) /A and our heuristic analysis suggests that the candidate rate 
function of Q\ (t, y) /A is given by 



sup 

•)6C 6 [0,T]x[0,oo) Jo 



0(t,y-t) 



d 2 



dt dy 



q(t, y)) dy-ip N lo 



3 e ^dF (y) 



o 



dt, 



it is then sensible to conjecture, making y = u + t, a representation based on g (i, it) = q(t, t + u) 
via 



sup 

6>(-,-)eCi,[0,T]x[0,oo) 



sup 

0( v )eCi,[O,T]x[O,oo) Jo 



6(t,u) 

e(t,u 



d 2 

dt dy 
d 2 



q(t, u + i) I du — i/>at I lo. 
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(u) 



dt 

e e(t ' u) dF {u 
(6) 



This representation, in turn, suggests that for the rate function to be finite at q (•), one might need 
to impose as a necessary condition the existence of d 2 q (t, u) jd 2 u. Nevertheless, as we shall see in 
our examples, one might have a finite- valued rate function even in cases in which dq (t, •) /du is not 
even continuous for every value of t £ (0, T). 



2.3 Construction of an Auxiliary Continuous Process 

In order to prove Theorem [l] we introduce an auxiliary approximating continuous process, Q\, 
which shall be shown to be exponentially equivalent to the process of interest Q\ in the uniform 
norm. The construction of Q\ will be based on polygonal interpolations, so it will be convenient to 
introduce some notation. 

First, given (t,y) and (t',y') where t / t' we write Q\ (t,y) Q\ (t',y') to denote the straight 
line that joins the points (t,y,Q\ (t,y)) and (t',y',Q\ (t',y')) in the associated three-dimensional 
space. 

Now, given a sample path of the process Q\(-), consider the set {ti,...,t m } of points cor- 
responding to either arrivals or departures in the interval [0,T] (in increasing order); and put 
to = and t m+ \ = T. First let us consider Q\(t,-) for a fixed time t G {to, ...,t m+ i}. Let 
{yi (t) , y n (t) (t)} be the set of discontinuities of the function Q\ (t, •) (recall again that Q\(t, •) is 
a right continuous non-increasing step function). Interpolate using straight lines forming the seg- 
ments Q x (t, 0) o Q x (t, yi (t)), Q x {t, yi (t)) O Q\(t, y 2 (t)), . . . , Q\(t, y n ( t )-i (*)) <-> Qx(t, y n{t ) (t)). 

The next step is to join the end points of these straight lines to the end points of adjacent 
(suitably matched in the time axis) end points of straight lines in order to form segments of adjacent 
planes. In order to do this matching note that for each successive ij and iffi, either Q\(ti + i, •) has 
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one less discontinuous point than Q\ (ti,-) (i.e. a departure occurs at ij+i) or one more discontinuity 
point (i.e. an arrival occurs at tj+i); the exception is the last segment from t m to i m +i = T, where 
there might be no difference between the number of discontinuity points between Q\(t m ,-) and 
Qx (t m +i, •)• Note that batch arrivals are not possible since the interarrival times are positive. 

According to the notation introduced earlier for discontinuity points, y\ (ti) , y n (tA (U) are 
the discontinuous points of Q\(U,-) with corresponding values Q\(U, y\ (ti)), Qx(ti,y2 (ti)), 
Qx(ti,y n (ti) (**))• We wm explain how to joint discontinuity points of Q\ (ti,-) with those from 

Qx (U+i, •)• 

Suppose a departure occurs at time U+i. Then we can label the discontinuous points of 
Q\(t i+1 , •) as yi (t i+1 ) , ■■■,y n ( ti+1 ) (U+i), with n (t i+1 ) = n(ti) - 1. We form a set of straight lines 
Qx(ti,0) Q\(t i+1 ,0) o Qx(U,yi (U)) o Q\(U+i,Vi (U+i)) Qx{U,V2 (*»)) «-> Q\(t i+1 ,y 2 (t i+1 )) o 
••• ^ <5A(ij+i,2/n(t l+1 ) (ti+i)) O Qx(ti,y n (u) (U)) in a zig-zag manner; together with another set of 
straight lines Q\(ti,0) -R- Qx(U,yi (U)) <->•...-<-> Qx(ti-,Vn(t l ) (U)), and also the set of straight lines 
Qa(*i+i,0) o QA(*i+l>yi ^ ■ • • -H- Qa(^+i, 2/n(ti +1 ) These three sets describe a series 

of adjacent triangular planar sections which jointly form a continuous surface. 

Similarly, suppose that an arrival occurs at time U+i. Then we can label the discontinuous 
points of Q\(U+i, ■) as Qx(U+i,yi (U+i)), . . ., Q\(U+i,y n (t i+1 ) (k+i)), with n(t i+1 ) = n (ti) + l. We 
then form the set of straight lines Q\(t i+ i,Q) O Qx(U,0) o Q x (ti + i,yi O Q\(U,yi (ti)) O 

Q\(U+i,y2 (U+i)) ••• ^ Qx(ti,y n (u)(ti)) Qx(ti+i,y n (u +1 )(U+i))- Again, together with a 
second set of straight lines Qa(^O) f-> Qx(U,yi (U)) o . . . o Qx(U,y n ( ti \ (U)), and a third set 
of straight lines, namely Q A (^+i,0) O Qa(^+1) 2/1 (*i+i)) ^ • • • Qx(U+l, y n (t i+1 )(ti+l)) ■ These 
three sets of straight lines, once again describe a series of adjacent triangular planar sections which 
jointly form a continuous surface. The last time interval from t m to T is dealt with similarly, with 
perhaps one less triangle formed if n (t m ) = n (T). 

The continuous function (Q* x (t,y) : < t < T, y > 0) is defined by concatenating all these 
adjacent triangular planar regions as one varies U and ij+i for i G {0, 1, m}, and setting Q* x (t, y) = 
for the region where y is beyond the boundary of the last triangular plane i.e. beyond the 
lines Qx(ti,y n (u)(U)) Qx(U+i, y n {t i+1 )(U+l)), i 6 {0, 1, ...,m}. It is immediate from the previous 
construction, and the fact that Q x (t, •) is non increasing, that QX(t, ■) is also non increasing for each 

te[0,T]. 

Then, we define our auxiliary process Qx (t, y) for y > t via 

Qx(t,y) = Ql(t,y-t). (7) 

In order to define Qx (t, y) for < y < t < T, first let Nx (•) be the continuous process obtained by 
the polygonal interpolation of N x (■), so that N x (0) = and N x (A k /X) = N x (A k /\) for all k > 1. 
Then, define 

Qx (t,y) = Qx (t A y,y) + N x (t) -N x (tAy), 
analogous to @. Observe that 

||Qa-0a||z><4, (8) 
where || • \ \x> represents the uniform norm over the set T>. 

3 Bounded Service Times 

In addition to the assumptions imposed in Section [2] here we also assume that P (V n G [0, M]) = 1 
for M G (0,oo). 
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We define T>m = {(t, u) : < t < T, < u < M + T} and let C + (T>m) be the space of functions 
(x (i, u) : (t, u) G T>m) such that x (•) is continuous in both components, x (0, u) = for u > 0, and 
x (t, •) vanishes at infinity. If in addition, x(-,-) is absolutely continuous, and dx(t,y) /dtdy = 
almost everywhere on < y < t < T, we say that x (•, •) G AC + (T>m)- 

Our initial goal is to obtain a large deviations principle for (Qx/X : A > 0) as A — > oo on the 
space (C+(T>m), IHId m ); we then will use Q to obtain the corresponding large deviations principle 
for (Qx/X : A > oo). 

We start by deriving a large deviations principle in the topology of pointwise convergence. The 
proof of this result will be given at the end of this section. 

Lemma 1. Let X consist of all the maps from T>m to M, and we equip X with the topology of 
pointwise convergence on T>m- Then Q\/\ satisfies a large deviations principle with good rate 
function I(q) defined by 



sup 

0(v)eC[O,T]x[O,Af] Jo 



T 



M+t / pp. \ / / pM 

0(t,y- t) I ? (t,y)J dy - ^ N (log (J e 9 ^dF (y) 



dt 



(9) 

if q(-) G AC+(T>m), and I(q) = oo otherwise. Here C[0,T] x [0, M] denotes the set of all continuous 
functions on [0,T] x [0,M]. 



In order to lift the large deviations principle indicated in Lemma [T] to the uniform topology we 
need the following result on exponential tightness; we shall also give the proof of this result at the 
end of this section. 

Lemma 2. Q\/\ is exponentially tight in C + (T>m) equipped with the topology of uniform conver- 
gence. 



Using the previous two lemmas we are ready to state and prove the main result of this section, 
which is a version of Theorem [T] for the case of bounded service times. 

Theorem 2. Q\/\ satisfies a large deviations principle with good rate function defined in ^ under 
the uniform topology on T>m- 

Proof. Since the domain of /(•) is a subset of C+ (T>m), and Qx/X G C + (T>m) with probability 
1, the large deviations principle in Lemma ft] holds in the space C + (T>m) with pointwise topology, 

(Lemma 4.1.5 (b) in [3]). Since by Lemma 



Qx/X is exponentially tight in (c + (T>m) , IHIx> M ) the 



same large deviations principle holds in [C+ (T>m) , ll"llx> A /J (Corollary 4.2.6 in [3j) and the result 
follows. □ 



As a corollary of the previous theorem we obtain that (Qx/X : A > 0) satisfies a large deviations 
principle on (L +i0O (P M ), IHIz> M ) 

Corollary 1. The process (Qx/X : A > 0) satisfies a large deviations principle on (L +)00 (T>m), 
IHIx) ) with rate function I (•) defined in ([£]). 
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Proof. First we verify that Q\/\ and Q\/\ are exponentially equivalent according to Definition 
4.2.10 in [3]). Since the laws of (Q\/X, Q\/X) are induced by a separable stochastic process and the 
underlying topology is induced by the uniform norm, the set 

{oj:\\Q x /X-Q x /X\\v m >v} 

is Borel measurable (see Remark b) following Definition 4.2.10 in [3]). Now recall that by the 
construction of Q\ that \\Q\ — Qx\\ < 4 a.s. Hence for any r\ > 0, 

P(\\Qx/X-Qx/X\\v M >v) = 

for large enough A. Hence 

limswp - log P(\\Q\/ \ - Qx/\\\v M > V) = 

A-s>oo 

The result then follows by applying Theorem 4.2.13 in [3]. □ 

3.1 Proofs of Technical Results 

Finally, we provide the proof of Lemmas [T] and [2} 

We start with Lemma [T] which takes advantage of the Dawson-Gartner projective limit theorem 
and thus requires that we obtain an auxiliary large deviations principle for finite dimensional objects 
defined via 

N x (ti) 

A u ( A ) = J2 J (%'-i < Ak / X + Vk ^ Vi) ( 10 ) 

fc=AT A (i l _ 1 ) + l 

= Q\(U,yj-i) - Qx(U,yj) - Qx(U-i,yj-i) + Q\(U-i,yj), 

for U-i < U, and yj-i < yj. 

Lemma 3. For = to < ti < h < ■■■ < t m < T and = yo < y\ < ... < y n < y n +i = T + M, 
(Ajj(A)/A : I < i < m, I < j < n + 1) possesses a large deviations principle with a good rate 
function 



sup 

9ij:l<i<m, 



P 6i 3 6i i ~J2 ^ N log Yl e&ijp (yi-i -u<V 1 <y j -u) \ du. 

l<i<n+l i=1 j=1 i=x J ti-x \ j=1 J 



Proof of Lemma\^ We use that i/'ivQ is continuously differentiable over M. Since Ui are non-lattice, 
the key renewal theorem implies that for any set of < to < h < t2 < • • • < t m , 

lim -\ogEex V \Y^9i{Nx{ti) - Nx{U-i)) \ = J2^n(0)(U - U-i) 

U=i J i=i 

for any 6 £ R. Then, from [1], the Gartner-Ellis limit of (Ajj(A) : 1 < i < m, 1 < j < n + 1) satisfies 

^ I m n+1 I 

A(9) = hm - log^exp } ^ %Ay(A) \ 

{■ i i i J 

m u I n+1 \ 

= J2 log e6ijp (yj-i -u<V 1 <y j -u)\du 

„■ i J ti— i \ - i / 



1=1 Jti ~ l \ j = l 
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is finite for any := (0j j : 1 < i < m, 1 < j < n + 1) . Moreover, for any < u < ti, 



J>_ 
d9in 



n+l 



-V'tv log e eifc P(y fe -i - u < Vi < y k - u) 

*> V fc=i 



n+l 



ip' N log Y e^Piyk-! -u<Vi<y k -u) 



fc=i 



e e ^P(yj-i — u < Vi < yj — u 
J2 n ktle e ^P(y k -i-u<V 1 <y k 



< max{ \ip' N (m&x{0 ik ,k = l,...,n + 1})| , \^' N (mm{9 ik , k = l,...,n + 1})| } 

which is uniformly bounded over a neighborhood of 6ij and < u < ij, fixing all other fy^'s. 
Therefore, 



V-iv log e 6 ^ +h P( yj -i - u < Vi < yj - u) + ^ e dlk P{y k ~\ -u<V x <y k -u) 



n+l 



V>tv log ^ e 9ik P(y k ^ - it < Vi < y fc 



k=l 



is also uniformly bounded one the same region. By dominated convergence theorem, we have 




A(6) 



n+l 



ip' N log ^2 e 6ik P(y k -i - u<V 1 <y k -u) 



Moreover, it is dominated by 



e^Pjy^i - u < V 1 < yj - u) 
EfeLi e e ^P{y k - X -u<V l <y k -u) 



du 



(U - U-i) max{\ip' N (max{9 ik , k = 1, . . . ,n + 1})|, |V^(min{6» ifc , A; = 1, . . . ,n + 1})|} < oo 

for any given B £ Since A(-) is finite and differentiate everywhere on M d , by the Gartner- Ellis 
Theorem for the case V\ = R d ([3], p. 52, Ex 2.3.20 (g)), {Ajj(A)} possesses a rate function 

m n+l m „i i I n+l \ 

su p yz yz ~yz ^ i°g ^^Pivj-i - u < vi < - «) 

8jj:l<i<m,l<3<n+lHj4 i^l Jt i-i \ j=l J 

We argue that the rate function is good. By [3] p. 8, Lemma 1.2.18, it suffices to show that 
(Ajj(A) : 1 < i < m, 1 < j < n + 1) is exponentially tight. Denoting || • ||i as the Li-norm, we have 
by ChernofF's bound 

IimA->oo^ log P (|| £s*j{\)/\\\ x > a) < BEa^ log P(N X (T) > aX) < -6a + ^ N {6), 
for any 9 > 0. Sending a + oowe then obtain 



lim^oolimA^oo- logP (|| A ij (A)/A|| 1 > a) = -oo, 

thereby obtaining exponential tightness and the goodness of the underlying rate function as claimed. 

□ 
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Proof of Lemma^ We will use the Dawson-Gartner projective limit theorem. Consider a collection 
of points in the plane of the form k = ((£$, y~) : 1 < i < m, < j < n), such that := to < ti < 
t2 < ... < t m < T and := yo < D\ < ■■■ < Vn- Moreover, we assume that y\ = t\ if < / < 
min (m, n). Let K be the union of such collection of sets k. Further, let {p k }kgk be the projective 
system generated by K. We will proceed to obtain a large deviations principle for the projections 
{Qx(t,y) I A '■ {t,y) S k). However, we will do this by first obtaining a large deviations principle for 
quantities Ajj (A) /A and then the large deviations principle for the projections follows using the 
contraction principle as the (Q\(t,y)/\ : (t,y) G k) will be shown to be continuous functions. Set 
y n +i = oo, so that Qx(t,y n +i) = for every t £ [0, T]. It is important to note, given the structure 
of the partition k, that if 1 < i < m, 1 < j < n, and i > j, then Ay (A) = 0. Now, similar to the 
definition of Ay (A) we define, for 1 < i < m and 1 < j < n + 1, 

Ay (A) = Qx(ti,yj-i) - Q\(ti,yj) - Qx(U-i, Uj-i) + Q x {U-i, yj). 

Once again, observe that Q x (t,y n+ i) = 0, and also if i > j, for 1 < i < m, 1 < j < n, we have 
U-i > Uj and therefore 

A^ (A) = Q x (yj-i,yj-i) + N x (U) - N x - (Qx{yj,yj) + N x (U) - N x (%•)) 

- (Q\(vj-i, yj-i) + n x - iv A (%_!)) + (Q x (yj,yj) + n x - iv A (%•)) 
= o. 

Moreover, clearly we have for 1 < i < to, and 1 < j < n 

i n+1 
Z=l r=j+l 

so indeed we have that (Qx(ti,yj) : 1 < i < m, 1 < j < n + 1) can be recovered as a continuous 
function of the A/ r (A)'s. Since Ay (A) — Ay (A) < 16, we have that 

^ I m n+1 j ^ I m n+1 j 

hm - log E exp ^ ^ % Ay (A) I = Urn - log £ exp Yl 9 *i^*i ( A ) \ • 

Consequently, from Lemma [3j the rate function for the projections represented by k (these projec- 
tions are denoted by p K (<j)) can be written as 

m n m / n+1 \ 

= sup ^2 V] <9y<5y (k) - ^ / log V] e 9i: >P(yj-i - u < Vi < yj - u) du 

{8 tJ :l<i<m,l<j<n+l} j=1 J=1 i=1 ./t^j I J=i / 

To possess a finite I(p K (q)), the quantity (/c) := q^yj-x) - q(U,yj) - q{U-i : yj^i) + q(ti-\,yj) 
must satisfy that 

kj («) = (11) 

for i > j, and 1 < i < m, 1 < j < n + 1; otherwise, if 5y (k) ^ 0, the rate function can be 
made arbitrarily large by picking 9ij = c x sgn(5ij (k)) with arbitrarily large constant c > for 
1 < j < i < m, as 

^ nog^e^ P^-i - u < Vi < yj - u) j du 
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is independent of 0y's that have j < i. In the representation of the rate function I(p K (q)) we have 
also used the fact that q(ti,t/j) = J2i<i r >j with g(0, yj) = 0, so the relation from the 5ij (k)'s 

to the q(ti,yj) is a one-to-one, continuous function, so that the contraction principle (Theorem 4.2.1, 
[3]) is invoked for the above representation for I(p K (q)). We want to show that swp KeK I(p K (q)) is 
equal to ([9]), and hence conclude the proof by Dawson-Gartner Theorem (see Theorem 4.6.1, [3]). 
Clearly it suffices to concentrate on functions g such that q(t, y) = whenever t>T or y>t + M 



given that we are assuming service times bounded by M. Note that the constraint (11 ) implies that 



for any q, in order that 1(g) < oo, we must have absolute continuity throughout < y < t < T 
and, moreover, that 

dq(t,y)/dydt = 

almost everywhere onO<y<t<T (see [3j p. 189). We now focus on q(t,y) that is absolutely 
continuous on [0, T] x [0, T -+ M] and have dq(t, y) jdydt = almost everywhere on < y < t < T. 
Observe that 



<% ( K ) = q(ti,Vj-i) - q(U,yj) - g(ti_i,y 3 -_i) + q{U-i,yj) 
ft, rvj Q% 

q(t,y)dydt. 



k-i Jyj 



_ 2 dtdy 



Regarding 6>(-, •) as a step function with jumps at = t\ < ti < • ■ ■ < t m < T and < yo < y\ < 
• • • < Vn < Vn+i = T + M, and denote S (D) as the set of all step functions on a given domain D. 
We can write 



sup I(p K (q)) 



sup 

0{;-)eS[a,T]x[o,T+M] Uo Jt 



t+M 



0(t,y) 



dydt 



q(t, y) dydt 



£ Viv (log e 9 ^dF (y - t)\ dt 



(12) 



To show that sup K /(p K (g)) > I(q) where /(g) is as defined in Q, note first that the set of step 
functions S ([0, T] x [0, T + M]) is dense in C ([0, T] x [0, T + M]), the set of continuous functions 
equipped with the uniform metric. So for any continuous function #(•, •) G C ([0, T] x [0,T + M]), 
we can find a sequence #&(•,•) G 5([0,T] x [0, T + M]) with \\9k — &\\[o,t]x[o,t+m] ~^ 0- Note that 
since 6 is continuous, it is bounded and so 6k is also uniformly bounded i.e. |#fe(i, y)\ < C for all 
and some C > 0. Consider 



o 



^(i,y) 



c) 2 



dydt 



q(t, y) dydt 



T / rt+M 

ip N log / e e(t ' y) dF(y - t) dt 



o 



with 9 £ C ([0, T] x [0, T + M]). We want to show that this can be approximated by the counterpart 
in 6 k g S([0,T] x [0,T + M]). Note that 



T rt+M 
it 



(t,y) 



2 



dydt 



q(t,y) 



dydt < C 



o Jt 



t+M 



d 2 



dydt 



q(t,y) 



since q is absolutely continuous. By dominated convergence we have 



10 Jt 



•Xt,y) 



d 2 

dydt' 



(t, y) dydt 



o Jt 



t+M 



0(t,y) 



dydt 



dydt < oo 



(t,y))dydt (13) 
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Similarly, since, as mentioned earlier \6k(t,y)\ < C, by the bounded convergence theorem we have 



t+M 



9k(t,y) dF ( y _ f ) 



t+M 



and so by the continuity of ipN (log (■)) we get 

*t+M 



9 ^y)dF{y-t) \ ^4> N lo. 



for any t. Furthermore, the obvious inequality 



t+M 



z e ^dF(y-t) 



rM j-M r-M 

/ dF(y) < / e e ^dF(y) < e C / dF{y) 
Jo Jo Jo 



(14) 



yields 



t/j N lo. 



t+M 



^v)f(y-t)dy 



< sup |Vjv(OI 

£e[-C,C] 



Hence yet another application of dominated convergence gives 



fl / /-t+M \ pT / f-t+M \ 

y ^at f log y e ek{t ' y) dF(y-t))dt^ j tp N I log y e e(t ' y) dF{y - t)\ dt (15) 

Combining (13) and (15) and using the expression in (12), we conclude that sup K /(p K (g)) > 7(g) 
(note a shift of variable y in (19])). For the other direction, consider 



T rt+M 



Jt 



0(t,y) 



2 



dtdy 



q{t, y) dydt - ^ N log / e e(t ^dF(y) dt 



M 



now with 6 eS ([0, T] x [0, T + M]). Note that we can find a sequence 9 k £ C ([0, T] x [0, T + M]) 
such that — > 9 pointwise almost everywhere and that 9 k is uniformly bounded; this sequence 
can be found, for example, by convolving 9 with a sequence of mollifiers (i.e. smooth kernels with 
bandwidth that tends to zero as k — >■ oo). Exactly the same argument as above would then yield 
sup K I(p K (q)) < I(q). Now, let q G C+(T>m) and suppose that q is not absolutely continuous. That 
is, it is not of bounded total variation in the sense of [3] p. 189. Then, for every 7 > there exists 

t x (7) < ... < t m (7) and y (7) < ... < y n (7) such that YT=i E"=i 



57- 
1.1 



> 7, where 



$3 = ?(*• (7) , V3-1 (7)) - i(U (7) , yj (7)) - (7) , yj-i (7)) + q(U-i (7) , (7))- 



Now observe that 
sup I(p K (q)) 

K&K 



sup 

K&K 



m n m r f. I n+1 

££m* (*)-£ / 1 ^ hog j^pfo. 
i=i j=i i=i y j=i 



u < Vi < yj — n) du. 



<5?-, and obtain 



Following [3J p. 192, we can select 9{j = sgn (5j- I for the partition introduced earlier that defines 



KdK 1 1 

1=1 j=l 
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Since 7 > is arbitrary we conclude that 



sup I(p K {q)) = 00 

K&K 



as required. 



Proof of LemmalR We want to prove that for any rj, 



□ 



lim — log P 

A— ¥00 A 



Qa(0,0) 



A 



> 77 



-00 



and 



lim liniA^oo i log P ( w I ^,5 j > rj 



-00 



where w(Q\/\, S) is the modulus of continuity of Qx/X with order S defined by 



w[Q\/\,5)= sup \Q\(h,yi)/\-Q\(t 2 ,y2)/\\. 
v J |ti-t 2 |<a 
\y\-yi\<& 

Recall that \\Q\ — Q\\\v M < 4 a.s., and that Q\(t,y) = Q\(t,y — t) for y > t and Q\(t,y) 
Qxiyji y) + N\ (t) — N\ (y) for < y < t < T. Therefore, it suffices to show that for any rj > 0, 



lim — log P 

A— too A 



Qa(0,0) 



A 



> rj 



-00, 



(where || • || denotes the supremum norm over [0,T] x [0, M]), also 



limlimA^ooilogP (w(^,S) > // 



-00, 



and finally that 



(16) 



(17) 



limliniA^ ilogP ( sup (iV A (i 2 )/A- N x (h)/\) > rj 

<>->0 ^ \0<t 2 -ti<<5 , 



-OO. 



(18) 



By our assumption that the system is empty, (16) is obvious. Condition (18) will follow as a direct 



consequence of our analysis of (17). Now, to prove (17) consider 

/ 



[T/S\ [M/S\ 



'(■(f.')>^E E p 

7 7 m=0 ra=0 



sup |Qa(*i, - Q\(h, 2/2)1 > Ary 

, 0<ti-t 2 <5, ti6(m<S,(ro+l)<5] 
V l?/2-J/i|<5, yie(n<5,(n+l)(5] 



It is best to proceed our analysis by keeping in mind the pictorial representation that we shall 
describe. One can represent the arrival and status of each customer in a two-dimensional plane, 
with x-axis representing the arrival time and y-axis the service time at the time of arrival. Under 
this representation, Q\(t,y) is the number of points in the triangle formed by a vertical line and 
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(a) Representation of Q(t, y) (b) Areas of Aff(m, n) and Af£(m, n) 

Figure 1: Illustrations for Q\(t,y) 



a 45° line passing through (t, y) in its northwest direction. Figure la depicts the shape of this 
triangle. Consequently, we have 



P 



( 

sup \Qx(h,yi) - Qa(*2, 2/2)1 > A77 

. 0<h-t 2 <S, tie(rfi<5,(m+l)(5], 
\ 1 1/2— Vl\<8, j/ie(n5,(n+l)(5] 



where 



< P(J\ff(m,n,X) + A/ 2 5 (m,n,A) > 77A) 



A^(m, n, A) = Q\(m6, (n - 1)5) - QaM, (n + 2)5) 



is the number of customers present at time mS who have residual service time between (n — 1)5 
and (n + 2)6, and 

JV A ((m+l)5) 

A/|(m,n,A)= ^ J(V5 > (n - 1)6), 

i=N x (mS)+l 

is the number of arrivals between m5 and (m + 1)5 which bring service requirements larger than 
(n — 1)5. Figure lb depicts the areas under which the points are included in J\ff(m,n,\) and 
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A/"! (to, n, A). Fixing 5 > and 9>0we take the limit as A — > oo in the following display, obtaining 
log£e e 



1 W K^(^(m,n,A)+A^ 2 4 (m,n,A)) 

A 

•Ar A (m5) 



r /N x (md) 

= -log£exp<j#f I(mS + (n-l)8-A i /\<V i <mS + (n + 2)8-A i /\) 

N x ((m+1)6) \ N 

+ Yl m>(n-i)5))\ 

i=N x (m8)+l / ) 

= - log E exp < y \og{e e P(m5 + (n - 1)5 - u < V { < m5 + (n + 2)<5 - u) + 1 

- P(m<5 + (n - 1)<5 - u < V t < m8 + (n + 2)5 - u))dN x {u) + log(e e F((n - 1)5) 
+ F{(n - l)5))[N x ((m + 1)5) - N x (m5)} j 

rmS 

/ ijj N (log(e e P(m5 + (n - 1)<5 - w < V < m5 + (n + 2)5 - u) + 1 
■/ o 

— P(m5 + (n - 1)5 - u < Vi < m5 + (n + 2)5 - u)))du 
+ ip N {\og{e e F({n - 1)5) + F((n - 1)5)))5 

=: ips(0;m,n) 

Fix 9 > 0. We argue that ^(0, m, n) — > as 5 — > uniformly over m, n. Indeed, for any m, n, 

T.5 



/ ip N (log(e 9 P(m5 + {n - 1)5 - u <Vi < m5 + {n + 2)5 - u) + 1 
Jo 



- P(m5 + (n - 1)5 - u < Vi < m5 + (n + 2)5 - u)))du 

< / ip N {\og{e e P{{n - 1)5 + u < Vi < (n + 2)5 + u) + 1 - P((n - 1)5 + u < V < (n + 2)5 + u)))du 
Jo 

< M^ N {\og{e e a{5) + 1 - a{5))) 

where a{5) := sup^^M] P( x < < x + 35) = o(l) as 5 — > oo by our assumption that the 
distribution of Vi is continuous and the fact that a continuous function is uniformly continuous on 
a compact set. On the other hand, we also have 

^(log(e e F((n - 1)5) + F((n - 1)5)))5 < ^ N {9)5 

for any m, n. Combining, we get 

MO, ™, n) < Mi) N {\og(e e u{5) + 1 - a{5))) + ^ N {9)5 (19) 

Now fix m and n. By Chernoff 's inequality we get 

P(Aff(m,n,\) +Afi(m,n,\) > rj\) < e -vO*+M6,™,n)\+o(x) 

and so 

hm^-^oo — log P (N i (m, n, A) +M 2 (m,n,X) > r]\) 

< — r\Q + ipg(9, m, n) 

< - V 9 + Mip N (log(e 6 a(5) + 1 - a(5))) + ^ N {9)5 
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by (19). Hence 



which gives 



liniA^oo j log P [ w 



Ox 
A 



_ [T/S] [M/S] 

< liniA^oo-log Yl p (Ni(m,n,\)+Ni(m,n,\)>ri\) 

m=0 n=0 

rj9 + M^ N (\og(e e a{5) + 1 - a(S))) + ^ N (0)S 



< 



lim lim sup — log P ( it; 



A 



,5} > rj) < -r]6 



Since 6 can be arbitrarily large, we conclude (17). Finally, condition (18) follows from the analysis 
of Af$(m,l,\)/\. 

□ 



4 Unbounded Service Times 



In this section, we will extend our result to unbounded service times. The main intuition of the 
extension beyond the bounded case is to justify that we can ignore in certain sense the customers 
who arrive with very large service time. Let us first introduce a suitable truncation scheme. For 
any K > and q G AC + (V) define 



for t e [0, T] and y := u + t > 0. 



t i-K 



d 2 



J y—w 



dsdz 



drdw 



s=w,z=r+w 



(20) 



The present 
customers with initial 









The present 




customers with 




initialservice 




times K 







Figure 2: Illustration for (f>fc(q)(t, t + u) 
Since q is absolutely continuous, <j>K (q) (t, u + t) is well defined. Moreover, the region over which 



the integration in (20) is performed corresponds to the triangular area depicted in light color in 
Figure [2} This region corresponds to the customers that are present at time t, have residual residual 
service time greater than y, and whose initial service time is less than K, as illustrated in Figure [2j 
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Moreover, for a sample path Q\, define Qx,K as the two-parameter process derived from Q\ by 
ignoring the arrivals with service time greater than K (one way to imagine is that they leave the 
system immediately upon arrival). Therefore, Qx,k is a two-parameter queue length process corre- 
sponding to an infinite server system with i.i.d. interarrival times following the law U = Y2i=i Ui/\, 
where G is a geometric r.v. independent of the U^s such that P (G = n) = F{K) n ~ l F(K), n > 1. 
It is easy to check that the arrival process corresponding to Qx,Ki i- e - by ignoring the arrivals with 
initial service time larger than K, satisfies the conditions in Section 2.1. The service time then has 
the distribution function Fk{x) = F{x)/F{K) for x £ [0, if]. We denote (Vn ,n = 1,. . .) as the 
sequence of service times in this modified system. 



Now recall the continuous version of Q\, denoted by Q\ constructed in Section 2.3 Moreover, 
define Qx,k to be the continuous approximation to Qx,K constructed in exactly the same fashion. 
In addition, for q £ AC + (T>k) define Ik (<?) as 

IT » - " {-WBy^ »>) * " {f{K)£ e9i '" )dFiy) 

(21) 

and set Ik (q) = oo otherwise, where is the infinitesimal moment generating function corre- 
sponding to the truncated arrival process. 

Theorem [i] yields that Q\,k/^i satisfies a full large deviations principle with good rate function 



sup 

6(t,-)eC[0,T]x[0,T+K] Jo 



Ik{-)- For q G AC+ (V) we shall also evaluate Ik (q) according to the expression (21). 

Since the geometric r.v. G is independent of the Ui's, we can compute the associated logarithmic 
moment generating function of the modified interarrival times 



K W(0) := k(6) + log 



1 - F(K)e K W 



and from which we solve that the associated infinitesimal logarithmic moment generating function 
of the arrival process is 

:= 4>N(log(F(K)e e + F(K))). 



dt, 



Plugging in the above expressions into (21), we have the following expression of Ik{q) 

sup f T \ [ T+K 6(t,y-t)(- J^-q(t, y))dy- tj; N (log (f{K) + f e e ^dF(y)) ) ] dt. 

6{;-)eC[o,T]x[o,T+K]Jo Ut \ at ay J \ \ J J J 

(22) 

At this point our strategy involves two steps. First, we want to show that Qx^k/^ and Q\,k/^ 
are exponentially good approximations as K /*• oo to both Q\/\ and Q\/X respectively. The second 
step consists in using this fact, together with the properties of Ik(q) as K f - oo and also properties 
of I (q) to conclude the identification of the rate function of Qx/X. 

So, to execute the first step we first define 

Nx(t) 

3=1 

that is, (t) is the number of arrivals with service time larger than K in the A-scaled system. 

Then we obtain the following result, which is proved at the end of this section. 
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Lemma 4. For any e > 0, 



lim lim^ - log P(ivf } (T) > Ae) = -oo. (23) 

K—>oo A 



Consequently, Qx,k/X and Qx,k/^ ar ^ exponentially good approximations as K /* oo to both Qx/X 
and Qx/X respectively. 

Using the previous lemma we obtain the following result. The proof is straightforward, but 
following our convention we shall give it at the end of the section. 

Lemma 5. The family {Qx/X : A > 0) satisfies a weak large deviations principle on C+ (T>) with 
rate function 

I* (q) := suplim^ inf Ik (z) ■ 

8>o {z--\\z-q\\T><°i 

We now extend the weak large deviations principle into a full large deviations principle with a 
good rate function using exponential tightness. 

Lemma 6. The family (Qx/X : X > 0) is exponentially tight on C+ (T>) and therefore it satisfies a 
full large deviations principle with good rate function /*(•). 

We proceed to show the identification I* (q) = I(q). We now collect useful properties that we 
will need to show this identification. 

Lemma 7. 

i) For any q such that I(q) < oo, we have I(4>K(q)) = Ik(4>k(q)) = Ik(q) /* I(q) as K — > oo; 
the notation Ik(q) I(q) implies that (Ixiq) '■ K > 0) is non decreasing in K and convergent to 

/(</)• ' ' ^ , 

ii) For any q such that I(q) = oo, and each M > 0, there exists a projection p K (following the 
notation introduced in the proof of Lemma [7]) such that, for large enough K, 

Ik(jpM) > M. 

Hi) Finally, with k from ii) there exists e > such that if q £ C+ (T>) and \\q — q\\x> < £ then 

hdpM) > M. 

We now are ready to prove the following important result of this section. 

Theorem 3. Qx/X satisfies a large deviations principle with good rate function defined in |i|) under 
the uniform topology on [0,T] x [0,oo). 

Proof of Theorem^ Given Lemma [6] all we need to show is that I* (q) = I (q). Suppose that q is 
such that I (q) = oo. Then, parts ii) and hi) in Lemma [7] imply in particular that for every M, there 
exists K, a projection k, and e > such that Ik {Pk (q)) > M for any ||<? — <?||x> < £• Consequently, 
we conclude, by using the monotonicity of Ik (q) as a function of K and taking subsequences, that 

J* (g) = suplim^ inf I K (z) = sup sup inf I K (z) = oo. 

5>0 {z--\\z-q\\ 7} <8} <5>0 _ftf>0 {^:||^-g|li3<o} 
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If I (q) < oo, then note that 



I* (q) = sup sup inf Ik (z) = sup sup inf Ik (z) . 

5>0 K>0{ z: \\ z -q\\-d^ S } K>0 <5>0 {« : l|2-9llp<<5} 

Further, observe that 

sup inf Ik (z) = sup inf _ Ik (4>k (z)) , 

s>o {^\\z-g\\ v <s} s>o {<t>K(z)-\\<t>K{z)-<t>K{q)\\ T ><&} 

Since Ik (•) is a rate function (in particular Ik (•) is lower semicontinuous) we have that 
sup inf I K ((f>K (z)) = Ik {4>k (<?)) 

<5>0 {<t>K{z)-\\<t>K{z)-<t>K{q)\\T> K <S} 

and then by part i) of Lemma [7] we conclude that sup^- >0 Ik {4>k {q)) = I (q), thus concluding that 
/* (q) = I (q) as claimed. □ 

We finish this section with the proof of Theorem [T] 

Proof of Theorem^ All we need to show is that Q\/\ and Q\/X are exponentially equivalent. This 
follows exactly as in the proof of Corollary [l] since \\Q\ — Qx\\v — 4. The measurability issue again 
is dealt with using separability The result then follows by applying Theorem 4.2.13 in [3]. □ 



4.1 Proofs of Technical Results 

We now provide the proof of the pending technical results. 

Proof of Lemma^ This result is a direct application of part a) in Theorem 4.2.16 in [H]. □ 

Proof of Lemma^ This is similar to the case with bounded service time, but the conditions for 
tightness are slightly different given that our domain T> is not compact. We must show that for any 
r], 7 > 0, we can choose small enough p > 0, such that for 5 < p, 

- log P(w(Q x /\, 6)>rj)< -7 

when A is large; this part is indeed basically the same as the case T>m- In addition, however, we 
also must show that for all rj > and every a > there exists K > such that 

P sup sup Qx (t, y) /A > ?? < exp (-Aa) . (24) 

\t£[0,T]y>K I 



Note that 



P(w(Qx/X,5) > v ) 

<P (w(Qx/\, 6) > v , \\Qx/X - Qx,k/M\ < |) + P (||Qa/A - Qx,k/M\ > |a 
<P (w(Qx,k/X, S) > |) + P (ivf > (T) > Xr}/2) , (25) 
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P I sup sup Q x (t,y) /A > 7? I < P ( ivf ) (T) > A77 



and 



s te[o,T] y>i^ 

By Lemma [4j for every 7 > can choose K large enough such that 

1 

A 



logP ( n[ K) (T) > Xt]/2 ) < -2 7 . 



for all A large enough. So, condition (24) is enforced and the second term in the sum in (25) is also 
appropriately controlled. Now, by a similar argument as in Lemma [2] in the previous section, for a 
chosen K, we have, for all small enough 5, 



A 



for large enough A. In sum, we get 



log P(w(Q x , K /\, 5) > I) < -2 7 . 



log P(w(Q x /\, 5) > 77) < -7 



for large enough A. Therefore, exponential tightness follows. It follows immediately that a weak 
large deviations principle and exponential tightness implies a full large deviations principle. The 
goodness of the rate function then is a consequence of exponential tightness together with the weak 
large deviations principle; see Lemma 1.2.18, p. 8, part b) of [3]. □ 

Proof of Lemma^ We start with part i), assuming that I(q) < 00. Since 



ff 2 



dtdy 



4>K(q)(t,y) 



& 2 



dydt 



4>K(q)(t,y) 



d 2 

dtdy 



<K+t}, 



we have immediately that I((/)k(q)) = Ik{<Pk{q)) = Ik(q)- It is obvious that Ik(q) is non-decreasing 
in K and that Ik (q) < I(q)- On the other hand, there exists 9 n G C^fO, T] x [0, 00) (the space of 
bounded and continuous functions on [0,T] x [0, 00)) such that 



I n (q) :-- 



T 



r(t, y -t) 



o 2 



dt dy 



q(t, y))dy-i} N [ lo 



9 n (t. 



y) dF{y) 



dt 



converges to I(q) as n — > 00. Since /(g) < 00, it follows easily that d 2 q(-)/dtdy is integrable over 
T>. Therefore, given that n {-) is bounded, 



JK 



O n (t,y-t) 



d 2 
dt dy 



q(t, y) dydt -> 



as K — > 00. Besides, for given n, as ipN is uniformly continuous on the bounded set [—M n ,M n 
where M n = sup \6 n {t,y)\ < 00, we have that 



N 



log (f(K) + j\° n ^dF( 



y) 



converges to 



dF(y) 
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uniformly on t G [0, T]. In sum, limx -5.00^(9) = I n (q) as IT — > 00. Therefore, there exists K n such 
that (g) > I n (q) — 1/n. Recall that (g) < Ik„(q) and consequently we obtain 



n 

Since Ik(,Q) increases in K, we have Ir(q) /* I{q) as claimed. 

For part ii), when I{q) = 00, there are two cases: a) q is not absolutely continuous, and b) q is 
absolutely continuous. Case b) in turn is divided into two subcases: b.l) d 2 q(t,y)/dt dy is not 
integrable over V, and b.2) d 2 q(t,y)/dt dy is integrable over V. We shall proceed to analyze all 
these cases now. For Case a). We can construct a projection p K with Ik(Pk(<i)) > M as we did in 
the proof of Lemma [T] For Case b), we have that q is absolutely continuous, but 



sup 

0(v)eC 6 [O,T]x[O,oo) Jo 



0(t,y) 



d 2 



dt dy 



q(t, y + t))dy-i/j N [ lo, 



y) dF(y) 



dt = 00, 



so we proceed to study case b.l): Assume that d 2 q(t,y)/dt dy is not integrable on T>. We shall 
assume that 

I {-dTdy^ y) ) dydt = °° (26) 

(if this integral is finite, then integral of the negative part must diverge and the analysis that follows 
next is identical). As in the proof of Lemma [TJ given a projection n induced by < t\ < ti < ... < 
t m <T and < y < yi < ... < y n +l, define 



5ij (k) := q(ti,Vj-i) ~ q(U,yj) ~ q(ti-l,Vj-i) + q(ti-l,Vj) 



(27) 



as long as yj-i > U-% (otherwise, if yj-i < U-x, Sij (/t) = 0). Then, from (26) and (27), it follows 
easily that for any M, there exists a partition k such that 



^8ij(K) >M + Tip N (l). 



Therefore, for large enough K, 



I K (p K (?)) > Yl 6i i ^ ~ T ^(l) > M. 



Now, for case b.2) suppose that d 2 q(t,y)/dt dy is integrable on V. We can find 0(t,y) such that 



3M < 



Jt 



6(t,y-t) 



d 2 



dt dy 



q(t, V)\ dy -ipN lo, 



dF(y) ) )dt < 00. 



Following the same line of reasoning as in the proof of part i) we can conclude that there exists 
K > such that Ixiq) > 2M. According to Dawson-Gartner Theorem, Ik{([) = sup Ik(Pk((i)) 
where the supremum is taken over all projections restricted to {t G [0,T],0 < y < t + K}. As a 
result, there exists some projection p K such that Ik(Pk(q)) > M and hence we are done. 
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Now we turn to part iii). So far we proved that for any q and M > 0, we can find a projection p K 
such that Ik(Pk(q)) > 2M. As discussed in the proof of Lemma [TJ we have 

m n+1 m „(. / n+1 \ 

SU P E E %%(«)-E / ' # lQ g E e 9ii P(Vj-i - u < < -u)\du> 2M. 

{9 ir .l<i<m,l<j<n+l} i=1 ~l i=1 Jti-! \ j=1 J 

where o~ij{n) is induced by the projection p K . By definition, there exists some 0y such that 

m n+1 m I n+1 \ 

E E - E / 1 #M lo S E e ^ P (%-l - u < V i K) ^ Vi -u))du> 3M/2. 

i=i j=i i=i \ j=i / 

For all e > and q E B £ (p), we have 

I M^) I <4e. 
Hence for e = M/(8^2 i j \0ij\) and all g E B £ (q), we have 

m n+1 m „j. / n+1 



I(p K (q)) > E E - E / # log ^ e^P( % _! - u < < % - - «) ok 

i=l i=l i=l Jt *-1 \ j=l ) 

m n+1 m I n+1 \ m n+1 

^ E E - E / ' ^n ] [ log E e ^ p (%-i - « < ^i (x) ^ yj - u ) ) du - 4 ^E E i^i 

i=l j=l i=l Jti - 1 V 3=1 / i=l i=l 

> M. 

Thus we conclude the result. □ 

Proof of Lemma [7} Let (T) be the total number of arrivals from time up to T with service 
time longer than K, under the A-scaled system. Then following [1] (or as in the proof of Lemma [2]) 
we have that 

lim - ] O gEe 0N >! C) m = T^ N {\og{e e F(K) + F{K))). 

A— >oo A 

Chernoff's bound yields 

P{n[ K) (T) > Xe) < exp{-#Ae + XT ^ N {\og{e e F{K) + F(K))) + o(A)}. 



Hence 



limx^jlog P X (N { X K \T) > Xe) 
<-6e + Tifj N (log(e F(K) + F(K))) 



Letting K — > oo gives 

lim hmx^ log P X (N { X K) (T) > Xe) < -9e 



Since 6 can be arbitrarily large, the result follows. 



□ 
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5 Examples 

This section is devoted to two examples that apply the large deviations principle that we have 
developed in the previous sections. The first example is on the most likely path to overflow in a 
loss queue, while the second example is on the ruin of a large life insurance portfolio that embeds 
an infinite server queue with service cost. 

Example 1. (Finite-horizon maximum of queue length process for M/G/oo) Consider an M/G/oo 
queue with Poisson arrivals with rate A. Suppose that the service times have a density / (•) with 
respect to the Lebesgue measure. The system initially starts empty. 

We want to find the optimal large deviations sample path to attain the event {maxo<t<T Q\(t, t)/X> 
x}, for fixed T and x, as A — > oo; this event corresponds precisely to the event of observing a loss 
in a queue with Ax servers, no waiting room, starting empty. Note that g (q) := maxo<«T q(t, t) is 
a continuous function under the uniform norm, so the contraction principle is directly applicable. 

We impose the condition that F(t)dt < x. This condition implies that the probability for 
the queue to reach Ax decreases exponentially fast as A — > oo (Such condition will be clear when 
we solve the constrained optimization in a moment). 

To proceed, let us first observe that iPn(6) = e 8 — 1. The maximization problem in ^ can be 
solved and the rate function is immediately recognized as 

which is easily seen to be a convex function of d 2 q(t,y)/dtdy. To find the optimal sample path 
amounts to solving the minimization problem 

Jo T JT (-£- y <l(t,y)) (log (-£- y q(t,y)/f(y - t)) - l) dydt + T 
subject to max < s <T Jq / s °° f-Jg^C*, y)j dydt > x, 

which is a convex optimization problem. The quantity Jq (—d 2 q(t,y)/dtdy) dydt is equal to 
q(s, s) when q is absolutely continuous, and q(s, s) represents the scaled queue length process at 
time s. 



To solve (28), we first consider a fixed s in the constraint and then optimize over s. When 



considering s fixed we replace the constraint in (28) by q(s,s) > x. Under this new constraint, it 



suffices to look at the time to s in the objective function, that is, we now solve 

&««,»)) (log 

subject to / " f n °° l—^q(t, y)) dydt > x. 



Jo ir v) ) Mog \-%^r\ ~ 1 1 *ydt + - (2g) 



The solution to (28) is then the optimal sample path from (29), among < u < T, that gives the 
smallest objective. 



We now consider (29). Introducing a Lagrange multiplier /x > 0, we minimize 
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By a formal application of Euler-Lagrange equations, we differentiate the integrand with respect to 
—d 2 q(t,y)/dtdy to get 

log( -^! ) y f^ )-^ = fory>u 
which gives 

d 2 _ u v_ f f(y-t) fort<y<u 
dtdy q[t!n 

> 1 (we replace by 

then implies 



dtdy \ nf(y — t) for y > it 

for some > 1 (we replace by another dummy /i for convenience). Complementary slackness 



{-d 2 q(t,y)/dtdy)dydt= / / l if{y-t)dydt = x, 
o -'m Jo Ju 

which in turn gives 

x 

^ = I, u F{t)dt 

(note that we have assumed j F{t)dt < x and so the condition > 1 is satisfied). As a result 

d 2 \ f(y~t) for t < y < u 

-m i(t ' V) = \f^i <ory>u. 

The optimal sample path q(t, y) leading to the constraint q(u, u) > x is given by 

* / ru roc x f( w _ g\ 

f(w - s)dw + / u - dw ) ds 
10 \JyAu JuVy J F(r)dr 

Transforming into q(t, y) = q(t, y + t) and some simple calculus gives the optimal sample path 



[ y+t - [ u - x fu-t F ( s ) ds 
q(t,y)= / F(s)ds- / F(s)ds^ " — , for y + t <u, 

Jy Ju-t Jn * V) (lr 

and 



iu-t Jo F(r)dr 

-y+t 



xH +t F(s)ds 
Jo F(r)dr 



In connection to our discussion about the direct rate function representation in terms of Q\/\ 
(see equation Q, in Section 2.2.1), one can check that dq(t,y)/dy is not continuous on the line 
y = t and therefore d 2 q(t,y)/dy 2 does not exists though I(q) is finite. 

Note also that the objective is 

f(y - t)dy + I" X / u ^~^ ( log ( ru - ^ - l) dy] dt + T 
nU ,U Ju J F(t)dt\ *\f u F(t)dtJ J y ) 

>( * )<ft+ ( l0g (^5)- 1 ) a - (30) 
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This is the rate function corresponding to the probability P(Q\(u, u) > Ax) = P(Q\(u, 0) > Ax), 
where Q\(u, 0) is the queue length at time u. This rate of decay is consistent with direct calculation 
using the fact that Q\(u,0) is a Poisson random variable with rate A F(t)dt, which gives 

P(Q(s) > Ax) = V e -^IoF(t)dt f x f U F(t)dt) jn\. 

n>Xx V J ° J 

For a consistency check, our result here can in fact recover the large deviations for the arrival 



process itself. If one changes the constraint in (29) to q(u, 0) = J* " J °° ( — ]^,Q(t, y)j dydt > x, the 



optimal value of ( |29| ) then becomes x[log(x/s) — 1] + s, which coincides with the exponential decay 
rate of P(Poisson(As) > Ax) as A — > oo. 

Figures [3] and [4] illustrate both the law of large numbers (i.e. the typical path) and the most 
likely path to the overflow event maxo< u <T Qx(u, u)/X > x for T = 1, x = 2. The underlying 
service time distribution is uniform in the interval [0,1]. We can see that the optimal path of 
Q(t,y) increases gradually over time to overflow at time 1. 
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(a) Q(t,y) for the most likely path to overflow (sur- (b) Q(t,y) for the unconditional path (surface) 

face) 

Figure 3: Surface plots of the asymptotic surface Q\(t,y) j A, as A increases, both an optimal (most 
likely) path leading to overflow, and the unconditional path. 



It is easy to see that since we assume Jq F(u)du < x, the rate function (30) is non-decreasing 
in s, and as a result an optimal time horizon is T. If the service time has bounded support [0, M] 
with M < T, then the selection of any time s £ [M, T] will give an optimal sample path. 

Example 2. (Insurance risk process) The net reserve of a life insurance company consists of 
the premium collected from policyholders, deducted by the benefit paid to policyholders in the 
event of deaths; often all these payments are discounted at zero in order to recognize the value 
of money in time. When policyholders arrive at the insurance company over time (an arrival is 
interpreted as the moment when a contract is signed), one can model the net assets of the insurer 
as a function of the underlying arrivals and death events of policyholders. Specifically, we shall 
assume that policyholders arrive according to a Poisson process with rate A, and that the time- 
until-death upon arrival of the policyholders are independent and identically distributed. Moreover, 
we assume that the time-until-death upon arrival has density /(•), distribution function F(-), and 
tail distribution F(-). The time-until-death in this setting can be thought as the service time in 
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(a) Q(t, y) for the most likely path to over-(b) Q(t,y) for the unconditional path 
flow (contour) (contour) 

Figure 4: Contour plots of the asymptotic surface Q\(t, y)/X, as A increases, both an optimal (most 
likely) path leading to overflow, and the unconditional path. 

the queueing context. We shall assume without loss of generality that the initial net reserve of the 
company is zero. 

It is often more convenient to work with the negative net reserve process, also known as the 
aggregate loss process, defined as the total benefit that the insurer has paid up to time t, minus 
the total premium received up to time t. For a policyholder who arrives at time Ai, and who dies 
at time A4 + Vi < t, the payoff by the insurer, discounted at time zero, is denoted h±(Ai,Ai + Vi); 
here Ai and Vi are the arrival time and time-until-death at the time of arrival of the policyholder. 
This quantity, h\{s,y), for y > s, captures the benefit paid at y minus the accumulated premium 
collected from time s to y. On the other hand, for a policyholder who has arrived prior to t, at 
time Ai, and who is still alive at time t, the payoff from the insurer to the policyholder is h 2 (Ai,t) 
(typically h 2 (Ai,t) will be negative as it represents premium that are paid to the insurer, so the 
payoff is negative). Here h 2 (s,t), for t > s, captures the premium accumulated from s up to the 
present time t, discounted to obtain the net present value at time zero. 

Consider, for instance, the setting of whole life insurance policies. That is, policies that pay a 
benefit b to the family of the policyholder, at the time of eventual death, in exchange of a premium 
which is paid at rate p continuously in time during all the time the policy was held, from arrival, 
up until the time of death. If the interest rate (or force of interest as it is known in the insurance 
setting) is constant equal to 8 > 0, then 

h x {s, y) = be' Sy - J\e~ Sr dr = be~ Sy - p{e~ 5s - e- 5y )/8, 

and ^ 

h 2 (s,t) = - [ pe- 5r dr = -p(e- 5s - e- &t )/8. 

J s 

The aggregate loss process, S\ (t), is represented as the net present value of the sum of the 
payoffs for all policyholders who arrive before t and it is given by 

N x (t) 

S x (t) =J2( J ( Ai + V ^ t ) MAi, ^ + Vi) + J(A + Vi > t)h 2 (A u t)) 

i=l 

ft ft ft f'OO 

= h 1 (s,y)dQx(ds,dy) + / / h 2 (s,y)dQx(ds,dy). 

JO Js JO Jt 
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We claim that S\ (•) is a continuous function of Q\ (•) under the uniform topology on T>. In order 
to see this, define D\ (t) to be the number of departures by time t, that is, 

D x (t) = N x (t) - Q x (t, t) = Q x (t, 0) - Q\ (t, t) . (31) 

Note that D\ (•) and N\ (■) are clearly continuous functions of Q\ (•). Moreover, we have that 

AW) rt rt 

i=1 Jo Jo 

and therefore 



V I(Ai + Vi< t) htiAi, A l + V i ) = / hi (s, u) D\ (du) N\ (ds) 

Jo Jo 



Sx(t)= h 1 (s,u)D x (du)N x (ds)+ / h 2 (s,y)Q x (ds,dy). 

Jo Jo Jo Jt 



Now, integration by parts shows that 



t rt 

hi (s, u) D x (du) N x (ds) 



o Jo 



t rt o2 







ds 

A similar development yields 



q 2 ft Q 

■ Q^-h l (s,u)Dx(u)Nx(s)dsdu- N x (t) J D x (u) —hi(t,u)du 
D x (t) f £-hi (s, t) Nx (s) ds + -Da (t) N x (t) h x (t, t) . (32) 



h2(s,y)Qx(ds,dy) = / Q\(s,y) — — dsdy - Q\(t,y) „ d V 



o Jt Jt jo 



dsdy J t dy 



+ f dk2 ^Qx(s, t)da - Qx(t, t)h 2 (t, t) . (33) 



It is now not difficult to see from (32) and (33) that indeed S\ (•) is a continuous function of Qx (■) 
in the uniform topology on [0,T] x [0, oo). 

Consider the finite-horizon ruin probability that the negative net asset of the insurer rises above 
the level Ax by time T. That is, the event {max tg r T i Sx (t) /A > x}. We wish to solve for the most 
likely path that leads to this event and therefore, applying our theory, we must solve the following 
convex calculus of variations problem. 

T roc (_J?_xU ..<\ ( ( 

subject to max < u <T Jo" (j? hi(t,y) (—^Q{t,y)) dy + h 2 (t,u) (^—^q(t,y)) dy\ dt > x 
Following the recipe of Example 1, we first consider 



*** So ir ( - v) Mog ^fS^r - 1 dydt + T 



ru roc 

mm 



fo U F (~& y ^ V)) ^g - 1 dydt + u 



subject to J " (f"hi(t,y) (-^g(*,3/)J dy + f™h 2 (t,u) (—^q(t,yj) dy^j dt > x 
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Introducing the Lagrange multiplier /j, > 0, we get 

_ JP_ -,. \ f" f(y - t)e^ hl ^ for t < y < u 

dtdy^ ' V > ~ \ f \y- t)e^^ for y > u 

When x is large, complementary slackness forces /x to satisfy 

lo (/ n ^^- t ) e ^ l( ^ ) ^( t '^ + ^- t ) e ^ 2(M)/l2 ^ u )) dt = x ( 34 ) 
for some \x > 0. Denote the integration on the left hand side by G(fi), then 

G'Ou) = QT /(y - ^^ft? (t, yjdy + F(u - t)e^^h 2 2 (t, u)\ dt > 0. 

Therefore, for given u, G(fj.) is monotone in \i. Besides, \G'(fi)\ — > oo as /j, — > oo. As a direct 
consequence, for any x large enough, equation (34) can be easily fit to many standard numerical 
solvers, and it admits a unique solution. Given fj>, the optimal sample path is given by 

q{t,y)= [ [ (--®^q{s,w)\dwds 



Jy 

t 



= 1(1 f{w-s)e^ u)hlis > w) dw+ I f(w-s)e^ u)h2{s ' u) dw)ds 

JO \JyAu JuWy J 

for y >t, and hence 

t I ru poo \ 

/ f(w- s)e^ u)hl ^ w) dw + f(w- s)e^ u)h ^ s ' u) dw ds 

\J(y+t)Au Ju\j(y+t) J 

for y > 0. Note that here we highlight the dependence of \i on u. Moreover, the rate function for 
the fixed-time probability is 

f( y _ t)e^ h ^\n{u)h 2 {t,u) - l)dy\ dt 

(35) 

The optimal time horizon u over < u < T is chosen to minimize (35). 

Now we consider a whole life insurance contract with benefit b = 1.5, continuous premium p = 1, 
zero interest rate and time-until-death which follows the uniform distribution on [0, 1]. Our goal is 
to compute the optimal sample path for ruin(x=10) before time T = 1. We solve the constraint 
equation (34) in Matlab and obtain the optimal u = 1 with /i = 2.251. 

In this case, we can compute the optimal path 

^ 1 ,// lilt — < i i i — i i i <,/,— <,J-i,t <,l*—it ^ _!;_:_:./ ^ 

and 



q(t, y) = 4(e^" ra - e^'^'^ - e^'^ + e^) + -e - ^"* - -j (e - ^ 1 - e - "), for y + t < 1, 



/x 2 /i /x 2 



q( t , y) = e-^-t-y) ( 1 — l^-m _ i-( e ^-^ - 1) ] , for y + t > 1. 

V m m / 

These optimal paths to ruin are shown in Figure [5} We just show the conditional paths, as the 
unconditional path are identical to the figures illustrated in Example 1. The optimal path here 
is qualitatively very different from that of Example 1. The value of Q(t,y) is the largest midway 
between time and 1. Intuitively, it is because it requires the smallest "energy", or distortion from 
the law of large numbers, at such time point in contributing to a large cash outflow from the insurer. 
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Time t 



(a) Q(t,y) for the most likely path to over-(b) Q(t,y) for the most likely path to overflow 
flow (surface) (contour) 

Figure 5: The surface and the corresponding contour plot of the asymptotic most likely path to 
ruin in a portfolio of life insurance policies. 
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